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■ Abstract. The equation of state of neutron star matter is examined in terms of the relativistic 
mean-field theory, including a scalar-isovector 8-meson effective field. The constants of the 

£vq , theory are determined numerically so that the empirically known characteristics of symmetric 

nuclear matter are reproduced at the saturation density. The thermodynamic characteristics of 
both asymmetric nucleonic matter and [3 -equilibrium hadron- electron npe-plasmas are studied. 
(*C) • Assuming that the transition to strange quark matter is an ordinary first- order phase transition 

CI described by Maxwell's rule, a detailed study is made of the variations in the parameters of the 

phase transition owing to the presence of a 8-meson field. The quark phase is described using an 
improved version of the bag model, in which interactions between quarks are accounted for in a 
one-gluon exchange approximation. The characteristics of the phase transition are determined 
for various values of the bag parameter within the range B £ [60,120] MeV/fm 3 and it is 
shown that including a 8-meson field leads to a reduction in the phase transition pressure Pq 
and in the concentrations um and uq at the phase transition point. 
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Qs . 1 Introduction 

O 

In addition to their independent, fundamental significance, studies of the structural characteristics 
and composition of the constituents of matter at extremely high densities and temperatures play 
an extremely important role in clarifying the physical nature of the internal structure and integral 
parameters of neutron stars. A quantum field approach in the framework of quantum hadrody- 
namics (QHD) provides a fairly adequate description of the properties of nuclear matter and of 
finite nuclei, treating them as a system of strongly interacting baryons and mesons. One theory of 
this type that has effective applications, is the relativistic mean-field theory [HO [3]. This theory 
yields satisfactory descriptions of the structure of finite nuclei [I] , the equation of state of nuclear 
matter [5], and the features of heavy ion scattering The parameters of the mean- field model 
characterizing the interaction of a nucleon with a, u, and p mesons can be self consistently deter- 
mined starting with empirical data on symmetric nuclear matter near the saturation density. This, 
in turn, leads to the possibility of obtaining equations of state for superdense, isospin-asymmetric 



1 Astrophysics, vol.52, No.l, pp. 132-150, 2009. Translated from Astrofizika, 52, No.l, pp. 147-164, 2009. 
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nuclear matter. In these studies it has been assumed that the masses of the scalar-isoscalar (a), 
vector- isoscalar (u), and vector-isovector {p) mesons and their coupling constants are independent 
of the density and of the values of the fields. In addition, the scalar-isovector <5-meson (ao(980)) is 
not included among the exchange mesons. 

Relativistic mean-field theory models have been constructed [8] assuming that the nucleon 
and exchange meson masses in nuclear media obey the Brown- Rho scaling law [9]. The results 
showed that including the density dependence of the mass leads to a more stiff equation of state 
for the matter. A scalar-isovector 5 -meson has been added to the scheme and a study of its role 
for asymmetric nuclear matter at low densities has been made in Refs. |10 | lil t f!2] . This approach 
has been used [131 CES] to study scattering processes for neutron-reach heavy ions at medium 
energies and the feasibility of forming a hadron-quark mixed phase during the collision process. 

This paper is a study of the equation of state for neutron star matter in terms of the relativistic 
mean-field theory and an examination of the variations in the parameters of the first order phase 
transition caused by including 5 -meson exchange. These results show how these variations will 
affect the integral characteristics and structure of hybrid neutron stars with quark matter cores. 

2 Lagrangian and thermodynamic characteristics of nucleonic 

systems 

The nonlinear Lagrangian density of an interacting multiparticle system consisting of nucleons 
and isoscalar-scalar u-mesons, isoscalar-vector w-mesons, isovector-scalar <5-mesons, and isovector- 
vector /9-mesons has the following form in qhd| 



C = ip N 



+ 2 



7 M (idp -g^u^x) - -S-pT^^Oc)) - (m N -g a <r(x) - g t l* lf ~T(x) y ) 
\ (d M <T(x)^a{x) - m a <r{x) 2 ^ - U{a{x)) + -m*t/ (x) - -9, ^{x)^ \x) 

\ (dj(x)d^(x) - m]-t(xf) + - Ut^lT (x) (1) 



where x = x^ = (t,x,y,z), cr(x), uj^{x), 5 (x), and ~ff >J \x) are the fields of the a, u, 5, and p 
exchange mesons, respectively, U(a) is the nonlinear part of the potential of the a-field, given by 

M 

U(a) = | m N {g a <jf + | (g a a)\ (2) 

/ i/j„ \ 

uijsf, m a , m^j, m$, m p are the masses of the free particles, ipN = is the isospin doublet for 

nucleonic bispinors, and f are the isospin 2x2 Pauli matrices. The symbol "— >" denotes vectors in 
isotopic spin space. This lagrangian, as in quantum electrodynamics, also includes antisymmetric 
tensors of the vector fields u^x) and given by 

( x ) = ( x ) ~ ( x ) ' $lfj,v (x) = d^p u (x) - d v pn (x) . (3) 



We shall use the natural system of units with Ti = c = 1. 
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Here g a , g w , gg, and g p in ([T]) denote the coupling constants of the nucleon with the corresponding 
meson. In the RMF theory, the meson fields a(x), u> p (x), 5 (x) and p P (x) are replaced by the 
(effective) fields a, Q^, 5, p^. Re-denoting the meson fields and coupling constants according to 

g a a = a , g w u) = u , g s 6^ = 5 , g p pf> = p , (4) 

(ga/m a ) 2 = a a , {g U} /m U} ) 2 = a u , (gs/m 5 ) 2 = a s , (g p /m p ) 2 = a p (5) 
and introducing the asymmetry parameter 

a = (n n - n p )/n , (6) 

the equations for the fields can be rewritten in the form 

a a (n sp (n, a) + n sn (n, a) - bm N a 2 - ccr 3 ) , (7) 

u = a^n , (8) 
5 = a 5 (n sp (n,a) - n sn (n,a)) , (9) 



a 



where 



p = --a p na, (10) 

k F (n)(l-a)V 3 

If m (a, 5) 2 

n sp (n, a) = — / p k dk , (11) 



k F (n)(l+a) 1/3 

y/k 2 + m* n (a, 5) 



n sn {n,a) = - 2 [ 7=^feS rx , k 2 dk , (12) 



7T 



.2^\ V3 



k F {n)=l — \ . (13) 

The effective masses of the proton and neutron are determined by the expressions 

m* (<r, 5) = ttin — a — 5, m* (a, 5) = tun — a + 5. (14) 

If the constants and a p are known, equations ([8]) and (fTUj) determine the functions uj (n) 
and p (n, a). Moreover, a knowledge of the other constants a a , as, b, and c makes it possible to solve 
the set of equations ([7]) , ([9]) , (fTT|) , (fT2|) in a self-consistent way and to determine the remaining two 
meson field functions a (n, a) and 5 (n, a). 

The energy density of the nuclear np matter as a function of the concentration n and the 
asymmetry parameter a has the form 

fcHnXl-a) 1 / 3 

\J k 2 + (rriN — a — 5) 2 k 2 dk 



1 

e{n,a) = 



7T 



+ ^2 / ^k 2 + (m N -o + 5) 2 k 2 dk + U{a) + -( — +— + *— + (15) 



7T- 



a 5 a p 
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where 

U (a) = - m N a 3 + ~ a 4 . (16) 
For the pressure of the nuclear matter we obtain 

k F (nXl-a) 1 / 3 

P(n,a) = -^ J (\/k F (n) 2 (1 - a) 2/3 + (m N - a - 6) 2 - y 'k 2 + (m N - a - 5) 2 ^ k 2 dk 
o 

+^2 y (J k F (nf (1 + a) 2/3 + (mjv - a + 5) 2 - ^ k 2 + (m N - a + 5) 2 ^j k 2 dk 



6(a) + l _£! + sd_£! + £ ! . (17) 

2 \ a CT a w a 5 a p I 



The chemical potentials of the proton and neutron are given by 
(J, p (n, a) = \J k F (n) 2 (1 - a) 2 ^ + (tun - cr - 5) 



2^ , 1 
+ OJ + -p, 



p n (n, a) = \fk~F {nf (1 + a) 2 ^ 3 + (m N - a + 5) 2 + uj - -p. (18) 



3 Determination of the constants for the model 

3.1 Empirical characteristics of saturated nuclear matter and the constants of 
the theory 

In order to determine the constants for the theory, a a , a w , as , a p , b, and c we can derive 
a system of equations relating these parameters to known empirical characteristics of symmetric 
nuclear matter at the saturation concentration no [17j . Given that the effective mass of a nucleon 
in symmetric nuclear matter (a = 0) at the saturation concentration no is related to the bare 
nucleon mass by 

m* N = ~/m N , (19) 

where 7 is a constant between 0.7 and 0.8, for the a field at the saturation concentration no we 
have 

(To = (1 - 7) m N . (20) 

Equations ([9]) and (|10p imply that 5q = and po = at the saturation concentration in 
saturated nuclear matter. Given the requirement that the energy e (n, a) jn per nucleon should 
have a minimum at n = n$ and a = 0, we obtain 

de (n, a) 



dn 



n = no no 
a = 



- m N + /o, (21) 
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where /q = B/A is the specific binding energy of the nucleus, neglecting the Coulomb interaction 
and finite size effects of the nucleus. 
Using Eq. ([TSJ, Eq. ^3TJ yields 

a UJ = — (m N + /o - J k F (n ) 2 + (m N - cr ) 2 ) . (22) 

The ojq field for symmetric matter at no, on the other hand, is given by 



^0 = a w «0 = m N + f - \J k F (n ) 2 + (m N - a ) 2 ■ (23) 
From the equation for the a field (J7|), we have 

k F (n ) 

CT o 2 f {m N -ao) 2 2 3 



> / /TTT7 r2 k dk ~ bm N^0 ~ ccr • ( 24 ) 

a CT 7T J ^/k z + (m N - a ) 2 

The energy density eq = no (tun + /o) f° r saturated nuclear matter at the saturation concen- 
tration no can be written in the form 



hp (no) 1 / 2 

:„ - ^ / V^ 2 + ( m ^v -<J ) 2 k 2 dk + -m N al + ^al + - ( ^- + nga w | . ( 2o) 



An important empirical characteristic which, in a certain way, couples the phenomenological 
constants of the theory, is the modulus of compression of the nuclear matter, which is defined as 



Substituting Eq. (fT5|) in Eq. (f26|) gives 



d 2 /s(n,a) 



■n 



n = n 
a = 



(26) 



K = 9a w n + 3- * '' 



k F (n ) + (m N - a ) 



n (m N - do) 2 1 



-9 u 2 JV 5 — — : . (27) 

k F (n ) 2 + (m N -a ) 2 1 2 F ^ o) fc4* _ 2 

[k z +(m N -ao) j 



In the semi-empirical Weizsacker formula the term for the specific energy of the asymmetry 
of the nucleonic system is given by 



--sym _ rp / \ 2 



n 

The coefficient of asymmetry energy, E S y m (n) , is defined as 

1 d 2 e (n, a) 



E S ym (n) a 2 . (28) 



Bsym (n) 



2n da 2 



(29) 

o=0 
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Using Eq. (|15p . we obtain the following expression for the symmetry energy at the saturation 
concentration of the nuclear matter, E^ m = E sym (no) 

^sym ~ a P ' 



N 6^ k F (n ) 2 + (m N - a ) ' 



1 M^v-^o) 2 1 (30) 

2 A* (no) 2 + ^-(70) 2 1 2 ("o) 

^7 + J 



«< ^ ^ q ( fc 2 + (m ^_ (To) 2)3/2 



3.2 Numerical determination of the constants for the theory 

In order to determine the constants for the theory we have used the following values of the known 
nuclear parameters at saturation: mjv = 938,93 MeV, 7 = m* N /mN = 0,78, saturation concen- 
tration of nuclear matter no = 0,153 fm -3 , specific binding energy /o = —16,3 MeV, modulus 
of compression K = 300 MeV, and E^y m = 32, 5 MeV. Equations ([20]) and ([23]) can be used to 
determine the <to and loq fields. Then Eqs. ([22]) . ([24"]) . ([25]) . ([27]) and ([30]) form a system of five 
equations for the six unknown constants, a a , a w , a$ , a p , b and c. It can be seen from Eq. ([30]) 
that including the interaction channel involving the isovector-scalar 5 -meson leads to a certain 
correlation between the values of as and a p . 



Table 1: Values of the Constant a p for Different Values of as 



as = (gs/m s ) 2 , fm 2 





0.5 


1 


1.5 


2 


2.5 


3 


a P = (9p/m p ) 2 , fm 2 


4.794 


6.569 


8.340 


10.104 


11.865 


13.621 


15.372 



Table 1 lists the values of a p for various values of the constant as- In order to clarify the 
role of the S -meson, in the following we shall take as = 2.5 fm 2 |llj . The absence of the 5 
interaction channel will correspond to an interaction constant as = 0. Note that the value used 
here, as = 2.5 fm 2 , is in good agreement with Ref. 18, where a microscopic Dirac-Bruckner- 
Hartree-Fock theory is applied to asymmetric nuclear matter and exotic nuclei in a study of the 
density dependence of the meson- nucleon coupling constants. According a plot of as as a function 
of concentration n in Fig. 2 of Ref. [18], the average value of as in the range n ~ 0.1 -f- 0.3 fm~ 3 
is on the order of 65 GeV -2 ~ 2.5 fm 2 . 

Table 2: Constants for the Theory without (crujp) and with (aujp5) a 5 -Meson Field 





a (j , fm 2 


a w , fm 2 


as, fm 2 


a p , fm 2 


b, fm 1 


c 


atop 


9.154 


4.828 





4.794 


1.654 • 10~ 2 


1.319 • 10~ 2 


atopS 


9.154 


4.828 


2.5 


13.621 


1.654 • 10~ 2 


1.319 • 10~ 2 



Table 2 lists the values of the parameters obtained by numerical solution of the system of five 
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Eqs. (i22j) . (f24l) . (j25j) . (|27f) and (|30j) without (crtup) and with (crwp^) the isovector-scalar 5 -meson 
interaction channel. 

4 Characteristics of {3 -equilibrium npe plasmas and the equation 
of state of neutron star matter in the nucleonic phase 

The values of the constants a a , a^, as, a p , b, and c for the relativistic mean-field theory obtained 
in the previous section (See Table 2.) can be used to calculate various characteristics of matter 
with an asymmetric proton - neutron composition (np -matter) and of (3 -equilibrium npe -matter. 
In terms of the relativistic mean-field theory, the lagrangian density for the npe -plasma is 

Cnm = C + ipe (i^d^ - m e ) ip e , (31) 

where C is the lagrangian of the system consisting of nucleons and au> p5 mesons (See Eq. ([TJ), ip e 
is the electron wave function, and m e is the electron mass. In this case, we find the energy density 
of the npe -plasma to be 

e NM (n, a, fjt e ) = e (n, a) + e e (u e ) , (32) 
where e (n, a) is the energy density of the npau p5 system defined by Eq. ([26]) . 

e e {pe) = — 2 J \jk 2 + m| k 2 dk (33) 
o 

is the contribution of the electrons to the energy density, and fj, e is the chemical potential of the 
electrons. For the pressure of the npe -plasma, we have 

1 / 2 2 \ 3 / 2 

P NM (n, a, fi e ) = P (n, a) + ^ [p e ~ m e) ~ e e (Me) • (34) 

Depending on the coefficient of surface tension, a s , it is known that a phase transition of 
nuclear matter into quark matter can occur in two ways [19] . It can either have the character of an 
ordinary first order phase transition with a discontinuous density change (Maxwell's rule) or mixed 
nucleon-quark matter can be formed with a continuous variation in the pressure and density |20j. 
In the second case, the condition of global electrical neutrality implies that, in order to determine 
the parameters of the phase transition and the equation of state of the mixed phase, it is necessary 
to know the equation of state of the (3 -equilibrium charged npe -plasma. To find the characteristics 
of the (3 -equilibrium, but not necessarily the neutral, npe -plasma, one has to solve the system of 
four equations (7)-(10) for specified values of the concentration n and asymmetry parameters a, 
and find the unknown mean meson fields a (n, a), u (n), 5 (n, a) p (n, a). Equations (18) can be 
used to determine the chemical potentials of the nucleons, p n (n, a) and u p (n, a), so that, using 
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the (5 -equilibrium condition, it is possible to find the electron chemical potential, and, ultimately, 
the energy density £nm and pressure Pnm of the (5 -equilibrium rape-plasma. 



He (n, a) = a n (n, a) - a p (n, a) 



(35) 



Figure 1 is a three dimensional plot of the energy per baryon, Eb(n,a) = enm/ii, as a 
function of the concentration n and asymmetry parameter a for the case of a j3 -equilibrium 
charged npe -plasma. The lines correspond to different fixed values of the charge per baryon, 
q = (rip - n e )/n = (1 - a)/2 - n e /n. 



E„. Mel 




2 a = (n -n )ln 




n, fm" 



Figure 1: Three dimensional representation of the 
energy per baryon Ef, as a function of the baryon 
number density n and the asymmetry parameter 
a in the case of a j3 -equilibrium charged npe 
-plasma. The upper surface corresponds to the 
"auipS" model, and the lower, to "crojp". The lines 
correspond to different values of the charge per 
baryon. 



Figure 2: The asymmetry parameter as a func- 
tion of the concentration n for a (3 -equilibrium, 
uncharged npe -plasma. The smooth curve corre- 
sponds to the "aujpS" model and the dashed curve, 
to "crwp". 



The thick line corresponds to (3 -equilibrium electrically neutral npe -matter. The lower surface 
corresponds to the "croup" model and the upper, to the "aujpd" model. Clearly, including a 5 -meson 
field increases the energy per nucleon, and the change is greater for larger values of the asymmetry 
parameter of the nuclear matter. For a fixed value of the specific charge, the asymmetry parameter 
falls off monotonically as the concentration is increased. Figure 2 is a plot of the asymmetry 
parameter as a function of the concentration n for the case of an electrically neutral " npe" -plasma 
for the "aujp" and "au>p5" models. It can be seen that including a 5 -meson field for a fixed 
concentration n will reduce the asymmetry parameter a. 

Figure 3 shows the effective masses of the protons and neutrons in a (3 -equilibrium uncharged 
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200 1 ■ ' ■ ' ' 1 

0.0 0.4 0.8 1.2 

n, fm" 3 

Figure 3: The effective nucleon masses as functions 
of the baryon concentration n for a j3 -equilibrium, 
uncharged npe -plasma for the "ctwjoJ" model. The 
dashed curve corresponds to the "aujp ,J model. 
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acop 




n - 


/" r— 1 


i . i 



3 kil 1 , 1 , L_ 

0.0 0.4 0.8 1.2 



n, fm" 

Figure 4: The concentrations of protons and neu- 
trons as functions of the baryon concentration n 
for a (3 -equilibrium, uncharged npe -plasma. The 
smooth curves correspond to the "puipS" model 
and the dotted curves, to the "aup" model. The 
dashed (straight) line corresponds to isospin sym- 
metric matter. 



npe -plasma as functions of the baryon concentration n for the "aujp5" model. Note that the 
effective masses of the protons and neutrons are the same in the "croup" model. Including the - 
meson mean field breaks the symmetry, in this sense, between the protons and the neutrons; the 
effective mass of the protons in this kind of medium is greater than that of the neutrons, i.e., there 
is a split in the values of the effective masses for the protons and neutrons. 

Figure 4 contains plots of the concentrations of the protons and neutrons as functions of the 
baryon concentration n for a /^-equilibrium uncharged npe -plasma. The dashed (straight) line 
corresponds to the case of isospinsymmetric matter. It is clear from this figure that the presence 
of a 5 -meson field reduces the neutron concentration and increases that of the protons. 

Our calculated equation of state for electrically neutral (3 -equilibrium npe -matter (the neutron 
star matter in nucleonic phase) using the "aajpS" model is shown in Fig. 5. Our equation of state 
(the segment of the curve labelled "MFT-aLupS") has been matched to the Baym-Bethe-Pethick 
(BBP) equation of state [21] in the region of normal nuclear densities. The Malone-Bethe-Johnson 
(MBJ) equation of state [22] is also shown for comparison. 
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Figure 5: The equation of state for neutron star 
matter in the nucleonic phase. The segment "MFT- 
aujpS" represents the results of this paper and 
"MB J," those of Ref. [22]. The region correspond- 
ing to nuclear - neutron (Aen) matter is described 
by the BBP equation of state [21] . 



5 The equation of state for a quark-electron ("iwise") plasma 

In order to describe the quark phase we have used an improved version of the MIT bag model 
[23], in which the interactions between the u, d, and s quarks inside the bag are accounted for in 
a one-gluon exchange approximation [21]. The quark phase consists of three quark flavors, u, d, 
and s and electrons that are in equilibrium with respect to weak interactions via the reactions 



d 



u + e + v e 



u + e 



d + v e 



u + e + v e 



u + e 



Since the v e and v e particles leave the system, the energy of the system decreases and the re- 
actions with neutrino emission continue until the condition [i v = holds for the chemical potential 
of the neutrinos. Then, the following conditions hold for the chemical potentials of the u, d, s, 
and e particles: 

Hd = Hs= H , Hit + He = H • (36) 

The following expression for the density of the thermodynamic potential Vlf of the quark flavor 
/ (/ = u,d,s) has been obtained in a framework of quantum hadrodynamics (QHD) [24J: 



Q f M = \ Hf^)-m) (/// - + \rn) In 



Hf + Ju f 



m 



-2— 

7T 



IMf + Juf 



jif 



m f 



2(/iJ-m? 
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- 3mj In 2 I — J + Qmj In ( — ) ( [ify/ftf - mj - mj In 



ml 



(37) 



where a s = <? 2 /47r, g is the QHD coupling constant, and p « m/3 « 313 MeV is the renormalization 
parameter. The quark concentrations are given by 



n f (a f ) 



H 2 — m 2 



f 



f J 'f 



1 2 n 

2 — mi — 2 — 



7T 



3m^ 



m 



The thermodynamic potential Q e and concentration of the electrons are given by 



(38) 



(At e ) = -^f / ~ \A 2 + k 2 dk , n e (/i e ) = — 3 JT^ ■ 



(39) 



The condition of electrical neutrality for a udse plasma is 

2 1 1 

-n u - -n d - -n s - n e = . 



(40) 



Using the functions n u (/x,/x e ), nd(fJ-), n s (fx), and n e (u e ) from Eqs. (38) and (40), this 
equation makes it possible to determine the function /j, e (u) and, ultimately, the functions 
^ u (a0 , tt d (u), Q s (//), and tt e (a) . 

The pressure of a "udse" plasma for a given chemical potential \i is given by 



Pqm (At) = - E *ii (M) - £ , 

i=u, d, s, e 



(41) 



where B is the "bag" constant, which characterizes the vacuum pressure and ensures confinement. 
The energy density £qm and baryon concentration nQM of a "udse" plasma are given by 



£qm(a0 = E (^i + AW) + 

2=U, d, 5, e 



(42) 



"QM(/i) = (i« + n (i + /i s )/3. (43) 

Equations (41)-(43) give the equation of state of a quark-electron ("udse") plasma in the 
parametric form eqm(P) and uqm(P)- As in the case of an npe -plasma, the baryon chemical 
potential for quark-gluon matter is given by 



Vqm(.P) = (P + e Q M(P))/n QM (P). 



(44) 
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6 Phase transition to quark matter at constant pressure 



The modern concept of the phase transition between nuclear matter and quark matter is based 
on a feature of this transition, first noted by Glendenning [20\ 117] , to the effect that there are 
two conserved quantities in this transition: baryon number and electrical charge. The requirement 
of global electrical neutrality then leads to the possible formation of a mixed phase, where the 
nuclear and quark matter are, separately, electrically charged, while overall electrical neutrality is 
ensured by electrons (leptons). In the case of a phase transition of this sort, the energy density e, 
baryon concentration n, and chemical potential [i e of the electrons, as well as the pressure P, vary 
continuously. The question of whether the formation of a mixed phase is energetically favorable, 
given the finite dimensions of the quark structures inside nuclear matter, the Coulomb interaction, 
and the surface energy, has been examined elsewhere [19j [25j I2S1 EZ]. It was shown there that 
the mixed phase is energetically favorable for small values of the surface tension between the 
quark matter and the nuclear matter. In this paper we assume that the transformation of nuclear 
matter into quark matter is an ordinary first order phase transition described by Maxwell's rule. 
A separate paper will deal with the changes in the characteristics of the phase transition with 
formation of a mixed phase [20j when the contribution of a 5 -meson field is included, as well as 
the influence of these changes on the integral and structural parameters of hybrid stars. In the 
case of an ordinary first order phase transition, it is assumed that both nuclear and quark matter 
are separately electrically neutral and that at some pressure Pq corresponding to the coexistence 
of the two phases, the baryon chemical potentials of the two phases are equal, i.e., 

Vnm (Po) = Hqm (Pa) ■ (45) 

Note that the chemical potential per baryon in nuclear matter is given by 

liNM = (/%> n p + fi n n n + nf M ">) /n , (46) 

and in the case of neutral, (3 -equilibrium nuclear matter (because of the conditions n p — rie = 
and [ip = [i n — fie NAI ^), coincides with the chemical potential /u n for a neutron given by Eq. (29). 
In the case of a neutral, (3 -equilibrium quark-gluon plasma, the relationship between the baryon 
chemical potential and the chemical potentials of a d quark (/^ = fi) and an electron ( fj,e ) has 
the form 

» QM = 3fi-^ M \ (47) 

7 Numerical computations 

Table 3 lists the calculated phase transition parameters for the "acop5+ MIT" model examined in 
this paper at constant pressure (Maxwell rule) for 12 different values of the "bag" parameter B. 
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Table 3: Parameters of a Maxwellian Phase Transition for Different Values of the "Bag" Constant B 



B 




n N 


n Q 


Po 


en 


£ Q 


(JNM) 


(QM) 
He 


A 


MeV/fm 3 


MeV 


fill -3 


fm~ 3 


MeV/fm 3 


MeV/fm 3 


MeV/fm 3 


MeV 


MeV 




60 


965.9 


0.1207 


0.2831 


2.11 


114.5 


271.4 


99.14 


9.205 


2.327 


65 


999.7 


0.1787 


0.3161 


7.22 


171.4 


308.8 


138.0 


8.350 


1.728 


69.3 


1032 


0.2241 


0.3504 


13.84 


217.5 


347.9 


166.0 


7.588 


1.504 


70 


1038 


0.2312 


0.3564 


15.10 


224.9 


354.9 


170.2 


7.464 


1.479 


75 


1079 


0.2810 


0.4027 


25.55 


277.6 


408.8 


198.1 


6.613 


1.349 


80 


1119 


0.3276 


0.4525 


37.95 


328.8 


468.6 


221.9 


5.842 


1.278 


85 


1158 


0.3704 


0.5036 


51.51 


377.5 


531.8 


242.1 


5.173 


1.240 


90 


1194 


0.4089 


0.5541 


65.54 


422.8 


596.2 


259.0 


4.605 


1.221 


95 


1227 


0.4435 


0.6029 


79.56 


464.7 


660.4 


273.1 


4.125 


1.213 


100 


1257 


0.4746 


0.6497 


93.30 


503.3 


723.5 


285.2 


3.717 


1.213 


110 


1309 


0.5281 


0.7369 


119.5 


572.0 


845.4 


304.5 


3.066 


1.223 


120 


1354 


0.5729 


0.8165 


143.9 


631.7 


961.4 


319.5 


2.568 


1.240 



The quark masses are taken to be m u = 5 MeV, = 7 MeV, and m s = 150MeV, and the 
strong interaction constant, to be a s = 0.5. In this table is the baryon chemical potential at the 
phase transition point, tin and tiq are the baryon concentrations of the nuclear and quark matter, 
respectively, at the transition point, en and eq are the energy densities, He and the 
chemical potentials of an electron in nuclear and quark matter, respectively, and Po is the phase 
transition pressure. 

It has been shown [28] that for a first order phase transition the density discontinuity parameter 

A = e q /(e n + P ) (48) 

plays a decisive in the stability of neutron stars with arbitrarily small cores made of matter from 
the second (denser)phase. Paraphrasing the conclusions of that paper |28j, in the case of a hadron- 
quark first order phase transition, we have the following conditions. If A < 3/2, then a neutron star 
with an arbitrarily small core of strange quark matter is stable. On the other hand, for A > 3/2, 
neutron stars with small quark cores are unstable. In the latter case, there is a nonzero minimum 
value of the radius of the quark core for a stable star. Accretion of matter to a neutron star when 
A > 3/2 will lead to a catastrophic (discontinuous) readjustment of the star, with formation of 
a star that has a quark core of finite size. This sort of catastrophic transition can also occur in 
the case of a rotating neutron star that is slowing down, when the pressure in the center rises 
and exceeds the threshold value, Po- The process of catastrophic readjustment with formation of 
a quark core of finite radius at the star's center will be accompanied by the release of a colossal 
amount of energy, comparable to the energy release during a supernova explosion. The last column 
of the table lists the discontinuity parameter A for the various values of the "bag" constant B. 
The above mentioned catastrophic readjustment of a neutron star (during accretion of matter to 
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Figure 6: The phase transition temperature Pq as a Figure 7: Baryon concentrations of nuclear matter 

function of the "bag" constant B. The smooth curve (tin) and strange quark matter (rig) at the point 

corresponds to the " auipS" model and the dashed of a maxwellian phase transition as functions of the 

curve, to "auip". "bag" constant B. Notation as in Fig. 6. 

its surface or as its rotation slows down) corresponds to the first three versions of the equation of 
state listed in Table 3, for which B < 69, 3 MeV/fm 3 . 

Figure 6 illustrates the dependence of the phase transition pressure Po on the value of the 
"bag" parameter B. It is clear that including a <5-interaction channel leads to a reduction in Pq. 
Similar plots of the baryon concentrations of the nuclear (njy) and quark (uq) phases at the phase 
transition point are shown in Fig. 7. Evidently, including a scalar-isovector effective S -meson field 
reduces the baryon concentrations for both phases at the phase transition point. Then the density 
discontinuity parameter increases. Figure 8 shows the equation of state for superdense matter with 
a maxwellian phase transition calculated in our "aujpS + MIT" model for five different values of 
the parameter B. 

8 Conclusion 

In this paper we have studied the equation of state of superdense nuclear matter in terms of the 
relativistic mean-field theory, including a scalar-isovector 5 -meson effective field. The values of the 
constants for the relativistic mean-field theory that we have found have enabled us to calculate 
the characteristics of asymmetric nuclear matter, as well as of (3 -equilibrium npe -plasmas. The 
dependences of the effective masses of protons and neutrons on the baryon concentration n for 
given values of the asymmetry parameter have been studied and it has been shown that in an 
asymmetric nucleonic medium the effective mass of a proton exceeds that of a neutron. 
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Figure 8: The equations of state for superdense 
matter with a maxwellian phase transition calcu- 
lated in the "acupS" model for five different values 
of the parameter B. 



The dependence of the asymmetry parameter a on the baryon concentration for j3 -equilibrium 
npe -plasmas has been studied for different values of the electrical charge per baryon and it was 
shown that including a 5 -field reduces a. 

Assuming that the phase transition between nuclear matter and strange quark matter is an 
ordinary first order phase transition obeying Maxwell's rule, we have made a detailed study of the 
effect of including a S -meson field on the parameters of the phase transition. We have determined 
the phase transition parameters for 12 different values of the bag parameter within the range 
B E [60 ; 120] MeV/fm 3 and shown that including a 5 -meson field leads to reduction in the phase 
transition pressure P$ and in the concentrations and uq for coexistence of the two phases. Here 
the density discontinuity parameter A increases. A bag parameter B ~ 69, 3 MeV/fm 3 corresponds 
to the critical value A cr = 3/2. When B < 69.3 Mev/fm 3 the density discontinuity parameter obeys 
A > A cr and neutron star configurations with infinitely small quark cores will be unstable. 

This analysis shows that a scalar-isovector 5 -meson field leads to more stiff equations of state 
for the nuclear matter owing to splitting of the effective proton and neutron masses, as well as to 
an increased asymmetry energy. It is known that a good source of information on the rigidity of 
an equation of state for dense matter is measurements of the mass of compact stars. The mass of 
the compact star in a binary system associated with the PSR pulsar 51516 + 02B was recently 
measured and found to be M = 2,08 ± 0, 19M [29]. The existence of neutron stars with such 
high masses argues for a more stiff equation of state than the equation which yields the standard 
value of M = 1,44 M©. 

Evidently, the above mentioned changes in the equation of state of superdense matter and in 
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the phase transition parameters will lead to corresponding changes in both the structure and the 
integral characteristics of hybrid stars with strange quark cores. A separate article will be devoted 
to a study of the configuration of neutron stars of this type, calculated by integrating the system 
of Tolman-Oppenheimer-Volkoff equations based on the equations of state obtained in this paper, 
with and without the inclusion of a scalar-isovector effective 5 -meson field. 
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